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We argue that the intrinsic glassy degrees of freedom in amorphous solids giving rise to the thermal 
conductivity plateau and the "boson peak" in the heat capacity at moderately low temperatures 
are directly connected to those motions giving rise to the two-level like excitations seen at still 
lower temperatures. These degrees of freedom can be thought of as strongly anharmonic transitions 
between the local minima of the glassy energy landscape that are accompanied by ripplon-like 
domain wall motions of the glassy mosaic structure predicted to occur at T g by the random first 
order transition theory. The energy spectrum of the vibrations of the mosaic depends on the glass 
transition temperature, the Debye frequency and the molecular length scale. The resulting spectrum 
reproduces the experimental low temperature Boson peak. The "non-universality" of the thermal 
conductivity plateau depends on ksTg/huiD and arises from calculable interactions with the phonons. 

PACS Numbers: 65.60.+a, 66.70.+f, 63.20.Ry, 66.35.+a 



A multitude of phenomena are observed in all low tem- 
perature glasses that can only be explained by the exis- 
tence of excitations not present in crystals. At very low 
temperatures, the universal scalings of heat capacity and 
thermal conductivity suggest these excitations cannot be 
simply ascribed to the local molecular motions specific 
to each substance . Yet quantitative analysis of the 
experiments has generally relied on purely phenomeno- 
logical models that invoke unrelated types of excitations 
in different temperature regimes. At liquid helium tem- 
peratures the excitations of glasses are well described as 
two- level systems ||-|| , but at somewhat higher temper- 
atures the more mysterious "boson peak" appears in the 
heat capacity and the phonon mean free path falls pre- 
cipitously, as if new scatterers became active, leading to 
a plateau in the thermal conductivity. Several different 
ideas about these additional excitations have been ex- 
plored, ranging from harmonic excitations of a disordered 
lattice H to anharmonic local modes of a "soft potential" 
with distributed parameters 

In contrast, here we suggest that like the two level sys- 
tems Q these quantum excitations in glasses are intrin- 
sic and essential relics of the non-equilibrium character 
of the glass after it is prepared. We show that quan- 
tized "domain wall motions" connected with the mosaic 
structure of glasses predicted by the random first order 
transition theory of glasses Q , explain quantitatively the 
boson peak and conductivity plateau. 

The random first order transition theory of the glass 
transition suggests there is a dynamical mosaic structure 
in classical supercooled liquids |l0[ . This mosaic struc- 
ture is directly manifested in the dynamical heterogene- 
ity observed in supercooled liquids using single molecule 
experiments |HJ and nonlinear relaxation experiments 
p2f. Below T g , the mosaic is spatially defined by the 



molecular motions that were not arrested at T g , and is 
thus only dynamically detectable. We have shown these 
same motions at liquid helium temperatures would be 
quantized and can be described by the two state system 
phenomenology. This identification predicts the density 
of states arising from these motions and accounts for the 
previously unexplained universality of thermal conduc- 
tivity in the liquid helium temperature regime. As the 
temperature is raised, however, the theory indicates these 
excitations should lose their strict two level character. 
The predicted temperature of the crossover to multilevel 
behavior coincides with the temperature regime where 
further excitable modes have been needed to explain ex- 
periments. In this paper, we show explicitly how the 
number of extra excitable modes of the mosaic accounts 
quantitatively for the data, thus removing the need to 
invoke additional mechanisms, although other contribu- 
tions may well be present to some extent. The excitations 
in the boson peak, like the two-level systems, active at 
lower temperatures, turn out to be collective motions of 
many particles (~ 200) encompassing rearrangements of 
a typical domain frozen-in at the glass transition. The 
multilevel behavior of these domains can be pictured as 
involving the concommitant excitation of the "ripplon" 
modes of the glassy mosaic along with the transition be- 
tween local minimum energy configurations of a mosaic 
cell. 

According to the random first order transition theory 
of supercooled liquids ||, metastable configurations be- 
gin to last locally longer than a few vibrations just below 
a dynamical transition temperature Ta- Below this tran- 
sition a glass forming liquid samples exponentially many 
metastable configurations much higher in energy than 
those of the corresponding crystal, but which are long 
lived. Glassy dynamics is described by the activated mo- 
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tions between minima. The liquid can be visualized as 
a mosaic of domains which locally resemble low free en- 
ergy configurations that are separated from each other by 
frustrated higher energy interfaces, or domain walls. The 
length scale of the mosaic and number density of these 
walls is determined by the competition between the en- 
ergy cost of a wall and the entropic advantage of using 
the large number of configurations. As temperature is 
lowered the configurational entropy per particle, s c , de- 
creases so that the cooperative length grows. This leads 
to a larger activation barrier that eventually gives relax- 
ation times exceeding the laboratory time scale. The size 
and rate of motion of these mosaic domains can be cal- 
culated without adjustable parameters once the liquid's 
heat capacity jump at the glass transition is known [flo|| . 
In the classical regime this theory leads to predictions of 
the average barrier and distributions of barriers in quan- 
titative agreement with experiment [fL3"| . According to 
this microscopic calculation, at the T g that corresponds 
with one hour relaxation time, the size of the cooperative 
region £ is expected to vary only a little from substance 
to substance. A mosaic cell allows for at least one alter- 
native kinetically accessible configuration at T g . These 
alternate states therefore have spectral density ~ 1/T g 
per domain of size £ 3 ||. This agrees with the mea- 
sured density of two-level systems ||. While activated 
at T g , below T gi motion between these alternative states 
remains possible, switching to coherent resonant tunnel- 
ing at a temperature proportional but smaller than the 
Debye temperature: ~ (a/^)T D /2-K M. 




FIG. 1. (a) Schematic of a tunneling mosaic cell is given, 
with doubled circles denoting atomic tunneling displacements. 
The boundary's location is variable, as illustrated (in an ex- 
agerated manner) by the blue and red lines. The domain wall 
distortion amounts to populating ripplon states on top of the 
structural transition energy e, as shown in (b). All transi- 
tions exemplified by solid lines involve tunneling between the 
intrinsic states and are coupled linearly to the lattice distor- 
tion and contribute the strongest to the phonon scattering. 
The "vertical" transitions, denoted by the dashed line, are 
coupled to the higher order strain; they contribute only to 
the Rayleigh scattering (which is too weak to account for the 
plateau 114]). 

At low temperatures the two-level system excitations 



involve tunneling of the mosaic cells typically containing 
N* ~ 200 atoms. The tunneling path involves stagewise 
motion of the wall separating the distinct alternative con- 
figurations through the cell untill a near resonant state 
is found. At higher temperatures, other final states are 
possible since the exact number and identity of the atoms 
that tunnel can vary (see Fig-ga). These new configura- 
tions typically will be like the near resonant level but 
will also move a few atoms at the boundary, i.e. at the 
interface to another domain. These fluctuations of the 
domain boundary shape can be visualized as domain wall 
surface modes ( "ripplons" ) . They cost a surface energy 
that varies with the size of the domain and have kinetic 
energy consistent with the mass of the moving domain 
wall. It is not surprising that the ripplon's frequencies 
turn out to be proportional to Wd, the basic quantum 
energy scale in the system just as does the crossover tem- 
perature. A detailed calculation of the ripplon spectrum 
would require a considerable knowledge of the topology 
and statistics of the mosaic. At each temperature be- 
low Ta the domain wall foam is an equilibrium struc- 
ture made up of nearly flat patches. According to the 
random first order transition theory, the effective sur- 
face tension depends on the curvature and vanishes at 
large radius of curvature as cr(r) oc r -1 / 2 . To approx- 
imate the spectrum we notice that the ripples of wave- 
length larger than the size of a patch will typically sense a 
roughly spherical surface of radius R — £(3/47r) 1 / 3 . The 
surface tension of the mosaic has been calculated from 
the classical microscopic theory and is given by <j(R) = 
l(k B T g /a 2 ) log((a/d L ) 2 /ire)(a/Ry/2 where d L /ais 
the universal Lindemann ratio. By itself such a ten- 
sion would collapse the mosaic but this tension should 
be compensated by stretching the frozen-in outside walls 
of neighboring mosaic cells. This compensating effect 
can be approximated by an isotropic positive pressure of 
a ghost (i.e. vanishing density) gas on the inside. Cal- 
culating the frequencies of the surface eigen-modes of a 
hollow sphere, subject only surface tension, is a classic 
problem of mathematical physics fl5f| . Accounting for 
the unusual r dependence of the surface tension <j(r) oc 
,.-1/2 nl0c ii nes the standard result for the frequencies is 
by a factor of 9/8. The eigenmodes have frequencies 
Luf = ^(a/fhwR 2 ){l — l)(l + 2), where fhw — {d L /a) 2 pa 
is the domain wall mass per unit area |8) and p is the 
mass density of the glass. The l-th. mode of a sphere is 
(21 + l)-fold degenerate. Using k B T g ~ pc 2 s a 3 (d L /a) 2 [§, 
one finds then w t ~ 1.34 uj D (a/C) 5/4 \/( l " !)(' + 2 )/4 ~ 
0.15 u>d \/ (I — !)(/ + 2)/4. Because of the universality of 
the (a/£) ratio 10 1, uii is a multiple of the Debye fre- 
quency. Due to the material's discreteness, there are 
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no harmonics of higher than (£^N*) [(R - a/2) /R] ~ 
9.. 10th order, a relatively large number, which makes the 
continuum approximation plausible. The lowest allowed 
ripplon mode is I = 2, and has a frequency of ~ ITHz for 
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silica. This is in agreement with the Boson peak in the in- 
elastic neutron scatering data |l6| ). 1 = 1 corresponds to 
a domain translation and is accounted for by the under- 
lying phonon momentum non-conserving tunneling tran- 
sition itself. The I = is a uniform domain dilation 
relevant to aging below T g . The existence of the domain 
wall vibrations lets us visualize the multilevel character of 
the tunneling centers as exhibited at temperatures above 
the TLS regime. A schematic of the resultant droplet 
quantized energy levels is shown in Fig.|l]b. 

The (classical) density of energy minima of a domain 
is n(e) — k ^ T e £ / fcBTg j8[, where the energy zero cor- 
responds to the (high-energy) configuration frozen-in at 
T g . The negative e's correspond to some of the very nu- 
merous but mostly unavailable lower lying energy states, 
accessed by tunneling. If we refer to the lower state of 
the resonant pair as the local ground state, the excitation 
energy density becomes n(e) = k * T e~l e l/ fcgTg (e > 0). 
The partition function of a domain with an excitation en- 
ergy e, possibly accompanied by populating vibrational 
states on top, is given by Z e = 1 + e~ l3t JJ, Zf l+1 , where 
Zi = 1/(1 — e~P hu)l } is the partition function of an ith 
order ripplon mode. Apart from the excitation e, we as- 
sume here each ripplon has a nearly harmonic spectrum. 
Note, our energy level scheme automatically insures that 
the atomic motions near the thermally inactive defects 
will only contribute to the regular lattice vibrations. The 
specific heat per domain, calculated from Z t and aver- 
aged over n(e), is shown in FigJ^. 
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FIG. 2. The heat capacity per domain, as follows from 
the derived TLS + ripplon density of states, divided by 
T 3 . This includes the Debye contribution. The thin line 
neglects phonon coupling and has T 9 /Td = 4. The ex- 
perimental curve for amorphous silica from UjV, originally 
given in J/gK 4 , is shown multiplied by h 3 pc 3 (6-K )(£/a) 3 /k%, 
where we used uj d = {c s /a)(6n 2 ) 1/3 , {{/a) 3 = 200 §, 
p = 2.2g/cm 3 ,c s = 4100m/sec and T D = Z42K §. Using 
the appropriate value of ksTg /Tiujd = 4.4 for a-SiC>2 would 
place the peak somewhat lower in temperature than observed. 

In order to estimate the scattering by the domain 
tunneling accompanied by the ripplons and thus the 
heat conductivity, we need to know the effective scat- 
tering density of states, the transition amplitudes and 



the coupling to the phonons. Any transition in the do- 
main accompanied by a change in its internal state is 
coupled to the gradient of the elastic field with energy 
g ~ pc 2 s J ds ■ d(r), where d(r) is the molecular displace- 
ment at the droplet edge due to the transition B . There- 
fore any transitions between groups marked with solid 
lines in Fig.|l|b are coupled to the phonons with the same 
strength as the underlying (TLS-like) transition, shown 
to be g ~ ^ksTgpcla 3 B, and no selection rules apply 
for the change in the ripplon quantum numbers because 
of the strong anharmonicity. We do not possess detailed 
information on the transition amplitudes, however they 
should be on the order of the transition frequencies them- 
selves, just as for those TLS's that are primarily respon- 
sible for the phonon absorption at the lower T whose 
transition amplitudes are also comparable to the total 
energy splitting. 

Since all transitions couple equally to the phonons, 
we can now calculate the density of the scattering 
states. If e > 0, the phonon absorbing transition oc- 
curs from the ground state. The corresponding total 
number of ways to absorb a phonon of energy Tim is 
= JcT den ( C ) £{„, m } 5 ( huJ ~ t £ + Dim n imhui m ]) = 

l/k B T g £ {nim} 0(w - £ ;m n lm un m )e-W"-Xi m »»««««), 
where we sum over all occupation numbers of the ripplons 
with quantum numbers l,m (m — —I.. I). The transition 
may also occur from the higher energy conformational 
state (e < in Fig.|l|b), and we have to compute p(ui) 
from these states too. To find p(u>), we compute the 
cumulative density of states Ne{w) = J_ rfe — e < 

E ) = Jo den ( e )J2{n lrn } 5 ( hu} ~ [J2lm n lrnhuim ~ e]). In 

a calculation parallel to obtaining Zi, this can be evalu- 
ated using an integral representation of the step function 
9 along with the use of steepest descent. The scatter- 
ing from excited states is proportional to p exc (uj,T) = 
dE f(E, T) dN E (Lu)/dE, where f(E, T) = 2/{e^ + l) 
gives the appropriate Boltzmann weights. 

Accurate calculation of the heat conductivity requires 
solving a kinetic equation for the phonons coupled with 
the multilevel systems, which would account for satura- 
tion effects etc. We utilize instead a single relaxation 
time approximation for each phonon frequency. The 
Fermi golden rule yields for the scattering rate of a 
phonon with Tilu ~ fc^T the relation tJ 1 ~ [p( w ) + 

p exc (u>,T)]. The heat conductivity then equals k — 
iEojWp( w )^ c s' The mean free path cannot be less 
than the phonon's wave-length A (which occurs at the 
Ioffe-Riegel condition). We account for multiple scatter- 
ing effects by putting l m f p = c s r w +A. At high T, the heat 
is not carried by "ballistic" phonons, but rather is trans- 
fered by a random walk from site to site, as originally 
anticipated by Einstein [H for homogeneous isotropic 
solids. 

The coupling of the multilevel excitations to phonons 
leads to significant frequency shifts and damping of the 
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resonant transitions. To compute these coupling effects, 
we replace the discrete summation over the different un- 
coupled harmonics J2i I dtu S(uj — uii) by integration over 
"lorentzian" profiles J2i I ^_^J)p; +7 2 , where y w is a 
(frequency dependent) friction coefficient and u>i(u>) is the 
renormalized ripplon frequency, which has been shifted 
due to the corresponding dispersion effects. The total 
relaxation rate of a transition involving more than one 
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FIG. 3. The predicted low T heat conductivity. The "no 
coupling" case neglects phonon coupling effects on the rip- 
plon spectrum. The (scaled) experimental data are taken 
from Q for a-Si (k B T g /huj D ~ 4.4) and [0 for polybuta- 
diene (ksTg/huio — 2.5). The empirical universal lower T 
ratio l m fp/l — 150 ||, used explicitly here to superimpose 
our results on the experiment, was predicted by the present 
theory earlier within a factor of order unity 

mode is thus the sum of the inverse life-times of the 
participating modes. This would be entirely correct for 
a frequency independent 7, but should be still an ade- 
quate approximation at the low T end of the plateau, 
where the absorption is mostly due to single ripplon 
mode processes. The mode decay rate 7^ is the phonon 
irradiation induced drag equal in the lowest order of 
perturbation to 7^ = j^(uj/c s ) 3 ~ §f k B T g (u/ui D ) 3 , 
where we do not distinguish between the longitudinal 
and transverse phonons. The corresponding (frequency 
dependent) ripplon frequency shift required by Kramers- 
Kronig relation gives a renormalized frequency uii(u>) = 



, where the princi- 



3 7 rri 1 1 \3 r^c du' (u)' /u>c) c 

pal value integral (numerically ~ 1) is a slow (but sign- 
changing!) function of w and the cut-off frequency lo c . 
lu c is of the order (but greater than) (a/£)w_D, since the 
phonons with wave-length shorter than £ cause an in- 
creasingly smaller effective gradient of the phonon field as 
sensed by a region of size £. Although our approximation 
f° r luj and 0Ji(ui) should break down in detail for lj > u> c , 
it makes little difference computationally, as the damping 
is already very intense at this point. Both the frequency 
shift and the damping broaden the absorption peak, the 
former being quantitatively more important. While the 
peaks' positions are determined by the quantum energy 



scale ujd, the spectral shifts are proportional to T g . The 
non-universality of the ksTg/hujD ratio, which varies 
over the range of 2 to 5, thus yields a non- universal po- 
sition of the plateau. This is consistent with the experi- 
mentally observed material dependence in this regime. In 
Fig.||, we show the resultant heat conductivities for differ- 
ent values of the ratio kBTg/Tiojo, using the specific cutoff 
lo c = 1.8(<i/£)wd. The non-monotonic behavior exhib- 
ited by this approximation in the plateau region is likely 
an artifact of neglecting thermal saturation in the kinetic 
theory. Assuming heat transport is dominated by ther- 
mal phonons suggests p exc (uj,T) -> p exc (aj,~ Uw/Jsb)- 
This leads to a temperature independent l m f p (uj) and a 
strictly horizontal plateau. 

The coupling effects on the heat capacity can be ob- 
tained by replacing each Zi with the formula for a 
damped oscillator. The resulting bump in c/T 3 , shown 
in Fig.|[ is therefore also non-universal depending on 
kBTg/hujo, being smaller for organic glasses, as seen in 
experiment. 
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